Renormalization of NN Interaction with Lorentz-invariant Chiral Two Pion Exchange. 
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I. INTRODUCTION 



At long distances the Nucleon-Nucleon (NN) interac- 
tion can be written as (see e.g. [l| for a review) 



V(r) = Vi,(r)+V2.(r) 



(1) 



where, Vi7r(r) and 'V2Tr{r) represent the One Pion Ex- 
change (OPE) and the Two Pion Exchange (TPE) con- 
tributions to the potential, respectively. Up to power cor- 
rections in the distance r, one has 'Vmr{r) = 0(6""'"''''). 
Such an expansion makes sense, since there is a clear 
scale separation at long distances. Actually, the omitted 
terms in Eq. ([T]) represent contributions whose ranges 
are shorter than l/(2m7r) 0.7fni. These include three 
and higher pion exchanges, correlated meson exchanges, 
etc 0- In momentum space, the expansion of Eq. ([T]) 
parallels an expansion on leading low momentum sin- 
gularities, rather than a naive low momentum expan- 
sion. The systematic and model independent determi- 
nation of those potentials was suggested several years 
ago and pursued by many others 0,[li,[il[ill 

(for some reviews emphasizing different viewpoints see 
e-g- [H, [H, [il, [H, [3 and references therein). However, 
in any scheme the potentials in Eq. ([1]) become singular 
at short distances, so one must truncate or renormal- 
ize the potential in a physically meaningful way in order 
to predict finite and unique phase shifts and deuteron 
properties. This has been a subject of much debate and 
controversy in recent times and we refer the interested 
reader to the literature for further detail s fiol . [itI . [isl . ITol . 

M M MMMM M M M M M 

In previous works by two of us (MPV and ERA) [25^, 
[2^, [28, 29] the renormalization of NN potentials was stud- 
ied using chiral potentials based on the heavy baryon 
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formahsm (HB-xPT) [8|, llO|. In the present paper, we 
extend those ideas to the case where Lorentz-invariant 
chiral potentials are used instead. In this case, the rel- 
ativistic framework of baryon chiral perturbation theory 
(RB-xPT) proposed by Becher and Leutwyler [sl, [s^ is 
employed in the construction of the two-pion exchange 
(TPE) component of the NN interaction. The remark- 
able difference between HB- and RB-TPE potentials lies 
on the long distance behavior [sl, [H, ISy, Is^l, due to 
the analytic structure of the ttN scattering amplitude 
in the low energy region where the Mandelstam vari- 
able t is close to 4to^, m^r being the pion mass. Ac- 
tually, one appealing feature of the RB-TPE potentials 
is that the long-distance two pion effects are correctly 
described, so that important contributions at the expo- 
nential level ~ g-2m^r properly re-summed, unlike its 
heavy baryon counterpart. In this work we are also inter- 
ested on its different short distance behavior, which plays 
an important role on the renormalizationprogram of the 
NN interaction developed in Refs. [H [lajp, We 
disregard, however, exphcit A's (see e.g. [38j) and other 
intermediate state contributions, assuming that those de- 
grees of freedom have been integrated out. 

As we have already mentioned, the calculation of scat- 
tering and bound state properties requires specifying the 
NN potential at short distances, which turns out to be 
highly singular for the Lorentz-invariant case. Actually, 
a crucial issue in the present context regards the number 
of necessary counter-terms required by the renormaliz- 
ability of the S-matrix. In the sin gle channel situation, 
the results found in Refs. [1^ HElHj HI] in configura- 
tion space can be summarized as follows. If the potential 
is a regular one, i.e., r^|[/(r)| < 00, there is freedom to 
choose between the regular and irregular solution of the 
corresponding Schrodinger equation. In the first case the 
scattering length is predicted, while in the second case 
the scattering length becomes an input of the calcula- 
tion. Singular potentials fulfill r^|[/(r)| 00 and do 
not allow this choice. For repulsive singular potentials 
[r^t/(r) — > 00] the scattering length is predicted while for 
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TABLE I: Sets of chiral coefficients considered in tfiis work. 





Source 


Ciyvjex ) 


ca^uev ) 




Set f 




-0.81 


-4.69 


3.40 


Set II 


NN \W] 


-0.76 


-5.08 


4.70 


Set III 


NN [42] 


-0.81 


-3.40 


3.40 


Set IV 


NN [20] 


-0.81 


-3.20 


5.40 


Set r? 


This work 


-0.81 


-3.80 


4.50 



The case 



the scat- 



'-U{r) 



9 IS 



attractive singular potentials 
tering length must be given, 
very special and, for g < —1/4, yields ultraviolet limit- 
cycles [13, H^, Hol • For coupled channels one must diago- 
nalize first the coupled channel potentials and apply the 
single channel rules to the outgoing eigenpotentials. 

In order to avoid any possible misunderstanding we 
hasten to emphasize that our use of the word relativis- 
tic is in a narrow sense; we are only disregarding a 
naive heavy baryon expansion of the virtual nucleon 
states in the calculation of the potential and hence tak- 
ing into account important anomalous thresholds singu- 
larities (3^ . Issj . This is not the same as providing a 
fully relativistic quantum field theoretical solution to the 
the two-body problem by say, solving a Bethe-Salpeter 
equation or any two body relativistic equation. This has 
always been a problem rooted in the non-perturbative 
divorce between crossing and unitarity in few body cal- 
culations for which the present paper has nothing to say. 

The paper is organized as follows. In Sec. [TTl we give an 
overview of our formalism already used in Ref. |25l . [26l . 
[28l . [29j . The key aspects on the derivation of the Lorentz- 
invariant TPE potential are briefly mentioned and the 
main differences with respect to the heavy baryon for- 
malism are highlighted in Sec. IIIII The deuteron bound 
state is discussed in Sec. IIVI Our predictions for phase 
shifts are displayed in Sec. El Finally, in Sec.|Vl]we draw 
our conclusions. 



II. FORMALISM 



Along the Hues of Ref. [25|, M IM ISi we solve the cou- 
pled channel Schrodinger equation in configuration space 
for the relative motion, which in compact notation reads 



u"(r) 



U(r) 



u(r) = fc^u(r) . 



(2) 



The coupled channel matrix reduced potential is defined 
as usual, U(r) = 2/i„pV(r), where finp = MpM^/ [Mp + 
Mn) is the reduced proton-neutron mass. For j > 0, 



U(r) can be written as 



33 



(3) 



U]]{r) 



In Eq. jg) P = diag(Zi(/i-hl),.. .,;Ar(;Ar-l-l)) is the angu- 
lar momentum, u(r) is the reduced matrix wave function 
and k the CM. momentum. In our case, = 1 for the 
spin singlet channel {I — j) and A^ = 3 for the spin triplet 
channel, with h — j — 1, h — j and ^3 = j + The poten- 
tials used in this paper were obtained in Refs. [sl. IssLIsgII . 
in coordinate space. We outline the main issues of this 
potential in Sec. IIIII 



A. Long distance behaviour 

At long distances, we assume the usual asymptotic nor- 
malization condition 



u(r) ^ h(-)(r) -h(+)(r)S, 



(4) 



with S the coupled channel unitary S-matrix. The corre- 
sponding outgoing and incoming free spherical waves are 
given by 



h(±)(r) = diagihf;{kr) 



(5) 



with h^{x) the reduced Hankel functions of order I, 



hi {x) = xH^^^^^{x) ( /ig {x) 



), and satisfy the 



free Schrodinger's equation for a free particle. 

The spin singlet state (s = 0) is an un-coupled state 



S\ 



(6) 



while the spin triplet state (s = 1) comprises one un- 
coupled I = j state 



(7) 



and the two channel coupled 1,1' = j ± 1 states for 
which we use Stapp-Ypsilantis-Metropolis (SYM or Nu- 
clear bar) [isj parameterization 



^J-lj-l ^3-1 j+1 



cos {2ej)e^'^^-^ isin {2ej)e'^^^-^+^^+^'> 
i sin (2ej)e''(^i-i+^l+i) cos {2e^)e^''^'^+^ 



In the present paper zero energy scattering parameters 
play an essential role, since they are often used (see 
below) as input parameters in the calculation of phase 
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TABLE 11: The number of independent parameters (counter-terms) for the Lorentz-invariant baryon expansion potential (RB) 
different orders of approximation of the heavy baryon expansion (HB) potential. The scattering lengths are in fm'+' +^ and 
are taken from Nijmll and Reid93 potentials [43| in Ref. [44|. We use the (SYM- nuclear bar) convention, Eq. (Illf) . The stars 
(*) mean that the behaviour is very dependent on the chosen set of chiral couplings. 



Wave 


a Nijmll (Reid93) 


OPE 


HB-NLO 


HB-NNLO 


HB-NNLO 


RB-TPE 










Set I, II & III 


SetIV 




'So 


-23.727(-23.735) 


Input 


Input 


Input 


Input 


Input 


'Po 


-2.468(-2.469) 


Input 


— 


Input 


— 


-(*) 


'Pi 


2.797(2.736) 


— 


— 


Input 


— 


— 


'Pi 


1.529(1.530) 


— 


Input 


Input 


Input 


Input 


'Si 


5.418(5.422) 


Input 


— 


Input 


Input 


Input 


'D^ 


6.505(6.453) 


— 


— 


Input 


Input 


— 


El 


1.647(1.645) 


— 


— 


Input 


Input 


— 




-1.389(-1.377) 


— 


Input 


Input 


Input 


Input 


'D, 


-7.405(-7.411) 


Input 


Input 


Input 


Input 


Input 


'P2 


-0.2844(-0.2892) 


Input 


Input 


Input 


Input 


Input 




-0.9763(-0.9698) 


— 


— 


Input 


— 


-(*) 


E'z 


1.609(1.600) 


— 


— 


Input 


— 


-(*) 


'Fs 


8.383(8.365) 


— 


— 


Input 


— 


Input 


'Fs 


2.703(2.686) 


— 


Input 


Input 


Input 


Input 


'D, 


-0.1449(-0.1770) 


Input 


— 


Input 


Input 


Input 


'Gs 


4.880(4.874) 


— 


— 


Input 


Input 


— 


Es 


-9.695(-9.683) 


— 


— 


Input 


Input 


— 




-3.229(-3.210) 


— 


Input 


Input 


Input 


Input 




-19.17(-19.14) 


Input 


Input 


Input 


Input 


Input 


'F, 


-0.01045(-0.01053) 


Input 


Input 


Input 


Input 


Input 


'H, 


-1.250(-1.240) 






Input 




-n 


Ei 


3.609(3.586) 






Input 




-n 


'H, 


28.61(28.57) 






Input 




Input 


'Hs 


6.128(6.082) 




Input 


Input 


Input 


Input 


'Gs 


-0.0090(-0.010) 


Input 




Input 


Input 


Input 


'h 


10.68(10.66) 






Input 


Input 




E5 


-31.34(-31.29) 






Input 


Input 





shifts. Due to unitarity of the S-matrix in the low energy 
limit, /e — > we have 



(S-1);,; =-2iar,;A:' 



I'+i+i 



(8) 



with ai'i the (hermitian) scattering length matrix. The 
threshold behaviour of the SYM phases is 







(9) 






(10) 






(11) 



B. Short distance behaviour 

The form of the wave functions at the origin is uniquely 
determined by the form of the potential at short distances 



(see e.g. [46', '4?'] for the case of one channel and [H, [2^, 
l28, 29] for coupled channels). For the Lorentz-invariant 
chiral NN potential , one has 

r' 

which is a relativistic Van der Waals type force (see 
e.g. [4^ for the electromagnetic case). Note that this 
short distance behaviour without an l/M expansion is 
at variance with the non-rclativistic and 1 /r^ in the 
standard Weinberg counting based on the HB chiral ex- 
pansion. In the latter, the expansion around the limit 
Af — > 00 is built in the formalism, leading to a different 
behavior at r — > 0. 

For a potential diverging at the origin as an inverse 
power law one has 

U(r) - ^ , (12) 



4 



with Cn a matrix of generalized Van der Waals coeffi- 
cients and n > 2. One diagonalizes the matrix C„ by a 
constant unitary transformation, G, yielding 



,n-2 



G 



(13) 



with Ri constants with length dimension. The plus sign 
corresponds to the case with a positive eigenvalue (at- 
tractive) and the minus sign to the case of a negative 
eigenvalue (repulsive). Then, at short distances one has 
the solutions 



"i,±(?') 



u(r) ^ G 



(14) 



where for the attractive and repulsive cases one has 



Tl/4 



sm 



— +Vi 



a. 



i/4 



exp 



Ri 
2 [ r 



(15) 



(16) 



respectively. Here, (pi are arbitrary short distance phases 
which in general depend on the energy. There are as 
many short distance phases as short distance attractive 
eigenpotentials. Orthogonality of the wave functions at 
the origin yield the relation 



N 

E 

i=l 



[Uk,: 



r=0 



A 

E 



(17) 



where A < N is the number of the short distance attrac- 
tive eigenpotentials. Details on the numerical implemen- 
tation of these short distance boundary conditions can 
be looked up in Refs. [H, [H, [H, [2|. 



C. Numerical parameters 

The Lorentz-invariant chiral TPE potential is specified 
by the pion weak decay constant Z^, the nucleon axial 
coupling constant gA, the nucleon mass M^r and the pion 
mass TOtt. In addition, at the level of approximation that 
we are working, it is enough to consider the low energy 
constants ci, C3 and C4 which characterize ttN scattering. 
The corresponding RB potential is specified by the same 
parameters at N'^LO in the HB chiral expansion. 

In our numerical calculations one takes /^r = 92.4MeV, 
= 138.03MeV, 2^„p = Mn = 2MpAf„/(Mp + M„) = 
938.918MeV, gA = 1-29 in the OPE piece to account for 
the Goldberger-Treimann discrepancy and gA — 1.26 in 
the TPE piece of the potential. The corresponding pion 



nucleon coupling constant takes then the value gjrNN = 
13.083, according to the Nijmegen phase shift analysis of 
NN scattering [4^ . The values of the coefficients ci , C3 
and C4 used along this paper can be looked up in Table [D 
There we list several Sets which have been proposed in 
the literature [13, H^j EH, S^l as well as the one which 
will be used in the present work based on our analysis of 
deuteron properties below. 

Renormalization requires fixing some low energy pa- 
rameters while removing the cut-off. We take the values 
from the high quality potentials [i^, [H^l as have been ob- 
tained in Ref. [4J] for the Nijmll and Reid93 versions. 
We will use the Nijm II values for definiteness. As men- 
tioned earlier, the number of independent parameters or 
counter-terms requires a study of the attractive/repulsive 
nature of the potential at short distances. The result of 
such an analysis for all channels considered in this work 
is summarized in table|TT]for the different parameter sets. 
We also list the scattering lengths in all partial waves as 
determined in our previous work p^ . 



III. LORENTZ-INVARIANT TWO-PION 
EXCHANGE 

A series of papers [H, [H, [H, [s^l was devoted to 
the construction of the TPE component of the NN in- 
teraction, based on the Lorentz-invariant formulation of 
baryon chiral perturbation theory proposed by Becher 
and Leutwyler [32, 33] in their study of the TriV sys- 
tem. These authors showed that it is possible to obtain 
a consistent power counting in a theory with a heavy 
particle without resorting to an integration of the heavy 
degrees of freedom and expansion around the limit of 
infinitely heavy baryon (HB-^PT). The latter results 
can be recovered from the Lorentz-invariant formulation 
through an expansion in l/rriN. However, this proce- 
dure destroys the correct analytic behavior of the ampli- 
tude near the low energy region close to < = 4m^. The 
underlying reason comes from the anomalous threshold 
of the triangle graph [32] (Fig[T]) right below threshold, 
t — Arri^ — m^/mj^. In the heavy baryon limit this singu- 
larity is ignored (as it collapses to the normal threshold) 
and any l/rriN expansion of the triangle loop around this 
region will fail to converge. Note that the same trian- 
gle integral also appears in the TPE potential, with two 
pseudo-vector vertices in one nucleon and a Weinberg- 
Tomozawa seagull term on the other. 

In order to illustrate the problem let us consider the 
spectral representation of the triangle graph. 



TT JAmi [t - t) 



(18) 



where 



Ini7(i') = -i ^ arctan n 9 — (19) 



t'-2ml 
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k-q 



p-k 



p-q 



FIG. 1: The triangle diagram, which cannot be reproduced by 
the usual heavy baryon expansion close tot — 4m^ . The solid, 
dashed, and wiggly lines represent, respectively, the nucleon, 
the pions, and an external scalar source. 



In HB-xPT the argument x = 2mMyJt' - 4m2/(t' - 
2m J) is assumed to be of order q^^ , yielding the expan- 
sion arctana: — 7r/2 — l/a; + l/'ix^ + • • • . The first two 
terms read 



ie., terms dependent on the variable z = p + p', where 
p and p' are the initial and final CM momentum of the 
NN system. The relativistic loop integrals, which incor- 
porates the dynamics of the TPE, also depend on this 
variable in a non-trivial way. However, phenomcnologi- 
cally one learns that such terms are not relevant at low 
energies and a Taylor expansion in z is usually consid- 
ered^. In this case the Fourier transform can be carried 
out in an easier way (see, for instance, Refs. [Sll, [53|). 
Generically, in any spin-isospin channel and up to the 
considered order in the RB-expansion the potentials may 
be written as a function of at most second order in the to- 
tal momentum operator. In this paper we take the local 
approximation on the radial part of the potentials and 
keep only up to linear terms in the operators. The re- 
maining non-localities are fairly small all over the range 
of interest, which somehow justifies its exclusion and con- 
siderably simplifies our calculations ^. 



IV. THE DEUTERON 



j{-q )\hb = 

1 dt' 1 



it' ~ 2ml) 



2T:m^A[q) H — — „ 9x ^(9) 



2 2mNy/t'^4ml 
, (20) 



where q = \q\, and L{q) and A{q) are the usual HE loop 
functions, 



^^^^ ^ y/^ml +q'^ ylm^+lf + q 



2m7r 



., , 1 q 

A[q) — — arctan 

2q 2m7, 



(21) 



However, it does not take into consideration the case 
I a; I < 1, when t' gets closer to 4m^. This region, where 
the naive heavy baryon expansion fails, is responsible for 
the long distance behavior of the triangle diagram, as can 
be seen by its representation in configuration space. 



r(r) 



dt' 



47r J4„2 



(27r)3 



1 ^-^q■r 



Im7(t') 



Im7(i') . 



(22) 



Therefore it is clear that, in order to have a good 
asymptotic description of r(r), one needs a decent rep- 
resentation for Im7(i') near t' = 4m^. This is only pos- 
sible if one takes the triangle anomalous threshold into 
account, which cannot be provided by current versions of 
the heavy baryon formalism. 

The potential in configuration space is obtained 
through a Fourier transform of the potential in momen- 
tum space. There one faces the problem of non-localities. 



In the pn CM system the deuteron wave function is 
1 



w{r) 



u{r)ap ■ cr„ 



(23) 



with the total spin s = 1 and nis = 0,±1 and dp and 
CT„ the Pauli matrices for the proton and the neutron 
respectively. The functions u{r) and w{r) are the reduced 
S- and D-wave components of the relative wave function 
respectively. They satisfy the coupled set of equations in 
the ^Si —3 Di channel 

— u"{r) -f (r)u(r) -f Ue^ {r)w{r) = — 7^u(r) , 



-w"{r) 4- UEi{r)u{r) 



w{r) 



(24) 



with (73 5j^(r), UEiir) and U3]j^(r) the corresponding ma- 
trix elements of the coupled channel potential. We solve 
Eq. p4)) together with the asymptotic condition at infin- 
ity 



u{r) 
w{r) 



Aoe 



1 + — 



7r (7?')^ 



(25) 



^ For instance, in the heavy baryon potential non-localities show 
up only at 0(g4) (N^LO). 

^ A rough way of estimating this in coordinate space is by acting 
with the operator V^/m^ on the local function which multiplies 
the variable and compare with the local which multiplies the 
variable z°. In the region between O.Olfm and 2fm these func- 
tions behave indeed as l/r"^ but the ratio between the non-local 
and local contribution is at the level of 1 — 5%. 
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r (fm) r (fm) 



FIG. 2: RB-TPE deuteron wave functions, u (left) and w (right), as a function of the distance (in fm) compared to the HB-TPE 
and Nijmegen II wave functions . The asymptotic normahzation u e"^"" has been adopted and the asymptotic D/S 
ratio is taken ri = 0.0256(4) in the TPE case (for OPE i] = 0.026333). We use the Set IV of chiral couplings and the Set rj (see 
table HI). 



TABLE III; Deuteron properties for the OPE and the HB- and RB-TPE potentials. We use the non-relativistic relation 
7 = ^2iinpB with B = 2.224575(9). The errors in the OPE case are estimated by changing the short distance cutoff in the 
range 0.1 - 0.2 fm. The errors quoted in the HB-TPE reflect the uncertainty in the non-potential parameters 7, rj only. The 
errors quoted in the RB-TPE are estimated as in the OPE case, but changing the cutoff in the range 0.3 - 0.4 fm. The entry 
Exp. stands for Experimental and/or recommended values and can be traced from Ref. '53j. 



Set 


7(fm-^) 


V 




r-d(fm) 


Qd(fm2) 


Pd 


{r-^)(fm-i) 


(r-=^)(fm-^) 


{r--'){im-^) 


OPE 


Input 


0.02634 


0.8681(1) 


1.9351(5) 


0.2762(1) 


7.88(1)% 


0.4861(10) 


0.434(3) 


00 


HB Set IV 


Input 


Input 


0.884(4) 


1.967(6) 


0.276(3) 


8(1)% 


0.447(5) 


0.284(8) 


0.276(13) 


RBE Set IV 


Input 


0.03198(3) 


0.8226(5) 


1.8526(10) 


0.3087(2) 


22.99(13) % 


0.5054(15) 


0.360(4) 


0.333(13) 


RBE Set ri 


Input 


0.02566(1) 


0.88426(2) 


1.96776(1) 


0.2749(1) 


5.59(1) % 


0.4438(3) 


0.2714(7) 


0.215(3) 


Nijmll 


0.231605 


0.02521 


0.8845(8) 


1.9675 


0.2707 


5.635% 


0.4502 


0.2868 


00 


Reid93 


0.231605 


0.02514 


0.8845(8) 


1.9686 


0.2703 


5.699% 


0.4515 


0.2924 


00 


Exp. 


0.231605 


0.0256(4) 


0.8846(9) 


1.971(6) 


0.2859(3) 


5.67(4)% 









where 7 = V MB is the deuteron wave number {B is the 
deuteron binding energy), Ag is the s-wave normalization 
factor determined from the condition 

/ dr [u{r)^ + w{rf] = 1 , (26) 

and the asymptotic D/S ratio parameter is defined by 
77 = Ad /As- As we have mentioned already the RB- 
TPE potential displays a relativistic Van der Waals 
singularity at the origin. Thus, the discussion on whether 
or not the deuteron parameters 7 and 77 can be fixed inde- 
pendently on the potential depends on the short distance 
behaviour of the eigenvalues of the c oup led channel po- 
tential matrix. As discussed in Ref. [2^ [2^, [2^, [2§| also 
for the bound state case the number of independent pa- 
rameters coincides with the number of negative (attrac- 
tive) eigenpotentials at short distances. In the RB-TPE 
potential we are using here [34,, 35, [sBl we have only one 
independent parameter (see table which we take to 
be 7 or, equivalently, the deuteron binding energy. With 



such a choice 77 becomes a prediction in contrast to the 
HB-TPE where t] is an input. The outgoing deuteron 
wave functions are depicted in Fig. [5] for the RB-TPE 
potential and compared to the HB-TPE one. 

Let us analyze in more detail the cut-off dependence 
of observables in the present RB-TPE potential. Given 
the fact that in the ^Si Di coupled channel we have 
one attractive and one repulsive eigenpotential at short 
distances we may borrow from the previous discussion on 
OPE [25! where we refer for further details. The practi- 
cal way of introducing in this case a short distance cut- 
off Tc which selects the regular solution at the origin is 
by appropriately choosing the auxiliary boundary condi- 
tion at the point r = rc among many possible choices 
compatible with self-adjointness The precise choice 
may provide smoother limits and hence better conver- 
gence properties in the pre-asymptotic region. Actually, 
as we show in Appendix \^ we can estimate the size of 
the finite cut-off corrections in deuteron observables and 
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FIG. 3: Cut-ofF dependence of the asymptotic D/S ratio of the 
deuteron wave function, 77, as a function of the short distance 
cut-off Tc (in fm) for the RB-TPE potential for the auxiliary 
boundary conditions u(r^ = and u'ijrc) = 0. We use the 
Set T) of chiral couplings and the Set r) (see table|l]| designed to 
reproduce the experimental value r] = 0.0256(4) in the limit 
rc ^ 0. The convergence towards the renormalized value is 
?7(rc)— 77(0) ~ exp[— 2/5(ii+/rc)'] up to oscillations (see main 
text). 

hence their convergence rate towards the corresponding 
renormahzed values. The result is that up to some oscil- 
lations the convergence towards the renormalized value is 
exponential as — > 0, i.e. the convergence towards the 
renormalized value is ri[rc) — 77(0) ~ exp[— 4/5(i?+/rc)2] 
for the RB-TPE potential ( ^ l/r"^ ) as compared to 
'7(^c) -'7(0) exp[-2(i?+/rc)^] for the OPE potential ( 
^ 1/r^). Here i?+ is a characteristic short distance scale 
of the corresponding repulsive eigenpotentials. The anal- 
ysis of Appendix [X] also shows that, generally, the more 
singular the potential the faster the convergence. 

For illustration purposes we show in Fig. [3] the calcu- 
lated asymptotic D/S ratio of the deuteron wave func- 
tion, 77, as a function of the short distance cut-off r,. (in 
fm) for the RB-TPE potential when the auxiliary bound- 
ary conditions u(rc) = and u'(rc) = are considered. 
Note the clear and coincident plateau below rc = 0.4fm 
for both boundary conditions following a relatively rapid 
variation above this region ^. As mentioned above, this 
is a typical feature of a coupled channel singular poten- 
tial with one attractive and one repulsive eigenpotential 
which extends to all other deuteron properties. Actually, 
the situation strongly resembles the previously studied 
OPE potential which has a softer 1 /r^ singularity at the 



^ Other auxiliary boundary conditions such as to(rc) = or 
w'{rc) = not shown in the figure make the oscillations de- 
duced in Appendix |X] more visible, but the plateau region takes 
place also below Tc = 0.4fm and yields an identical numerical 
result as with u{rc) = and u'ijc) = 0. 



origin [2^; the main difference is that for OPE stabil- 
ity takes place at a shorter scale, Tc = 0.2fm, than in 
RB-TPE potential. In the present context it is also help- 
ful to remind that a short distance cut-off radius and 
a sharp momentum cut-off are inverse proportional to 
each other, rc = 7r/(2A) [s^l (the numerical coefficient 
depends on the particular regularization) . Thus plotting 
observables as a function of rc enhances the finite cut-off 
changes while long plateaus could be observed instead as 
functions of A. Actually, halving the short distance cut- 
off corresponds to doubling the momentum space cut-off. 
For instance, in RB (Set-77) the range rc — 0.3 — 0.4fm 
corresponds to the range A = 780 — lOSOMeV where ob- 
servables such as ?7 change less than 0.01%. Pinning down 
this error bar would be harder if the sharp or other mo- 
mentum space cut-off was used. 

Our results for several deuteron properties are shown 
in table mil and compared to the corresponding HB-TPE 
considered in our previous work [H m m, ia]. Some 
remarks concerning the errors quoted in table IIIII are in 
order. We provide the largest source of error in the cal- 
culation which we are able to quantify. Since we aim 
at renormalized results, we stop whenever the change of 
the cut-off causes no significant variation within a pre- 
scribed accuracy, which we take to be about an order of 
magnitude higher than the typical experimental or rec- 
ommended value uncertainty. Thus, the cut-off range is 
not necessarily the same in all cases. In general terms, the 
more singular the potential at short distance, the faster 
the convergence of the result towards the renormalized 
limit (see e.g. Appendix [X}. Thus, the toughest case 
is OPE, where we only have 1/r^ singularity. Conver- 
gence in this case is the slowest, therefore shorter cut-offs 
Tc = 0.1 — 0.2fm are needed. 

The value we obtain for rj for the parameter Sets of ta- 
ble [T] is slightly different from the experimental one, mak- 
ing the comparison with the HB-TPE case [1^ , where ry 
was a free parameter (two attractive short distance eigen- 
potentials), a bit misleading. In order to obtain an accu- 
rate value of rj it was necessary to readjust the low en- 
ergy parameters C3 and C4 to the values C3 = — 3.8GeV~^ 
and C4 = 4.5GeV~^, indeed, very similar to the values 
proposed by other authors [13, H, El, IHI (see table [l|. 
Actually, once we have reproduced rj we see a general 
and slight improvement in accuracy when going from the 
HB-TPE (where 77 is a free parameter) of our previous 
work [2^ to the present RB-TPE calculation (where 77 is 
predicted). Basically, RB-TPE produces a sharp predic- 
tion for 77 (with eventually no errors), whereas HB-TPE 
does not predict ry, so its 1% experimental uncertainty 
propagates to other deuteron observables at about a sim- 
ilar 1% level, which is still comparable or larger than the 
error in the quoted experimental or recommended val- 
ues. The conclusion in Ref. [2^ was that agreement was 
partly achieved because of this fuzziness in the theoreti- 
cal predictions. Of course, one should not overstress the 
possible accuracy of the present results as regards the 
systematic errors; the main point of our calculation is to 
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FIG. 4: Relativistic Baryon expansion (RBE) np phase shifts in the ^So, '^Po and ^Pi channels as a function the short distance 
cut-off radius Vc for the fixed laboratory energies Tl ~ 10, 50, lOOMeV. 



provide the general picture when more complete asymp- 
totic TPE effects are correctly taken into account. 



V. PHASE SHIFTS 

We come to the calculation of the np phase shifts. In 
practice, this requires a careful wave by wave study of the 
renormalized limit. As can be seen in tablelTTl all coupled 
triplet channels have one attractive and one repulsive 
short distance eigenpotential. On the other hand, 
almost all singlet and uncoupled triplet channels develop 
an attractive singularity at short distances (some- 
times depending on the parameter values). The only ex- 
ceptions we found are the ^Pi and "^Pq channels, the latter 
depending on the precise values of the ci_3_4 constants of 
the chiral potential. This fact determines not only the 
number of counterterms, but also the convergence pat- 
tern towards the renormalized result. It reaches stability 
for cut-offs ranging in the region Tc = 0.3 — 0.5fm, de- 
pending on the particular partial wave and also on the 
energy (see the discussion in the previous section ITVl and 
Appendix |^ . Phase shifts in coupled channels with one 
repulsive singular component have been computed with 
either auxiliary boundary conditions uoj,i=j_i(rc) = 
and u'q J i^j_-f^{rc) = for zero energy states and sub- 
sequent orthogonalization of the finite energy states by 
using a complementary boundary condition as described 
in detail in Ref. [l^ for the OPE case. It is important 
to realize that even though renormalization requires in 
principle to pursue the mathematical limit Vc —^ 0, con- 
vergence is achieved in practice by length scales which are 
not unrealistically small and in fact are rather reasonable. 
This can be seen in Fig.|3]where, for illustration purposes, 
some selected low phases are depicted as a function of 
the short distance cut-off for fixed laboratory energy val- 
ues. As one generally expects smaller values of Vc are 
needed as the energy is increased. The approach towards 
the renormalized value for any partial wave depends on 
the attractive/repulsive character of the singular poten- 
tial at short distances. So, the ^Sq and "^Pq channels are 
purely attractive and hence the finite cut-off corrections 



are 0{rc^^) up to oscillations [40*1. On the other hand, 
the ^Pi channel provides a repulsive case and hence finite 

cut-off corrections are 0(6'^" ). We remind that the 
smallest de Broglie wavelength probed in NN interaction 
below pion production threshold is A ^ 0.5fm. Thus, the 
RB-potential and the present renormalization construc- 
tion also implement the desirable a priori requirement 
that short distance details are indeed irrelevant for the 
description of low energy properties. 

In Figs. [5] and [H] we present the np (SYM- nuclear bar) 
renormalized phase shifts for the total angular momen- 
tum j = 0,1,2,3,4,5 for Spin Singlet and Uncoupled 
Spin Triplet and Coupled Spin Triplet channels respec- 
tively. There we compare the relativistic baryon ex- 
pansion (RBE) and the heavy baryon expansion (HBE) 
as a function of the LAB energy compared to the Ni- 
jmegen partial wave analysis [43, HOl- For definiteness 
we use the chiral constants ci, C3 and C4 of Ref. Q 
(Set IV), which already provided a good description of 
deuteron properties after renormalization [26| at NNLO. 
This choice allows a more straightforward comparison to 
the N^LO calculation of Ref. [20] with finite cut-offs. We 
also compare with the Set 77 which takes to the same 
value of ci and the readjusted values C3 = — 3.8GeV~^ 
and C4 = 4.5GeV~"'^ based on our improved description 
of the deuteron in Sec. IIVI Unless otherwise stated, the 
needed low energy parameters for these figures are al- 
ways taken to be those of Ref. [i^l for the Nijmll poten- 
tial (see Table |TT]) . As can be clearly seen, the RB-TPE 
with this Set-77 not only improves deuteron properties 
but also the phase shifts all over. Again, one should not 
overemphasize this agreement, but it is rewarding to see 
that there is a general trend to stability and improve- 
ment in some channels (such as ^D2,^Pi, '^Po) when the 
RB-TPE potential is considered, while the quality of de- 
scription is not worsened in other channels. At the same 
time one should stress, however, that generally speaking 
this potential needs much less counter-terms as the cor- 
responding HB counterpart (about a half). Actually, in 
the Lorentz-invariant potential case one has at most one 
free parameter for channel as compared to the three pa- 
rameters for channels in the coupled triplets. This is due 
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FIG. 5: np (SYM-nuclear bar) Spin Singlet and Uncoupled Spin Triplet phase shifts for the total angular momentum j = 
0, 1, 2, 3, 4, 5 for relativistic baryon expansion (RBE) and for the heavy baryon expansion (HBE) as a function of as a function 
of the LAB energy compared to the Nijmegen partial wave analysis [431. ISOl]. 



to the attractive-repulsive short distance character of the 
coupled channel RB-TPE as compared to the attractive- 
attractive HB-TPE potential in these coupled channels. 
Again, we remind that the RB-TPE provides the correct 
analytic behavior of the exchange of two pions at large 
distances when A and other excitations are not explicitly 
considered. 



VI. CONCLUSIONS 



In the present paper we have analyzed the renormaliza- 
tion of all partial waves for NN scattering and the bound 
deuteron state for the Chiral Two Pion Exchange Poten- 
tial computed in a relativistic baryon Expansion. Our 
main motivation has been to consider a potential where 
the asymptotic TPE effects are consistently taken care of. 
This gives us some confidence that long distance physics 
is faithfully represented by a common exponential fall- 
off factor e"^™"'". At the level of calculation considered 
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FIG. 6: np (SYM-nuclear bar) Coupled Spin Triplet Phase shifts for the total angular momentum j = 0, 1, 2, 3, 4, 5 for relativistic 
baryon expansion (RBE) and for the h eavy baryon expansion (HBE) as a function of as a function of the LAB energy compared 
to the Nijmegen partial wave analysis \43i l50l | . 
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here A and other excitations are not explicitly taken into 
account; within a relativistic baryon expansion the con- 
tribution of this degree of freedom to the NN potential 
has never been computed. In addition, in the relativis- 
tic baryon expansion there appear non-local terms in the 
potential which involve the sum of initial and final mo- 
mentum space operators which fortunately turn out to be 
rather small all over the range as compared to the local 
contributions to the potential. This simplifies the analy- 
sis tremendously since standard coordinate space meth- 
ods can be used to solve the non-perturbative scattering 
problem and deuteron bound state and their subsequent 
and necessary renormalization. As we have repeatedly 
stressed along the paper this Lorentz-invariant potential 
presents a 1/r^ singularity at the origin which demands 
renormalization in order to get a finite and unique re- 
sult when the TPE potential is assumed to be valid all 
over the range from the origin to infinity. This can be 
done by introducing a number of (potential independent) 
counter-terms and consequently physical renormalization 
conditions must be specified. In practice they are fixed 
to the values of threshold parameters, mainly scattering 
lengths at zero energy. Actually, we have noted that the 
number of necessary counter-terms is drastically reduced 
when the Lorentz-invariant baryon potential is compared 
to the heavy baryon expanded TPE potential, while both 
potentials are specified by the same parameters. Thus, 
less input is needed to predict the NN phase-shifts. Al- 
though the precise number of counter-terms depends on 
the parameters of the potential we find that, typically, 
for the channels with total angular momentum j = to 
5 we need 13 in the Lorentz-invariant case as compared to 
about 27 in the HB potential. Actually, it is noteworthy 
that with about a half of the counter-terms the overall 
agreement is improved. This is particularly striking in 
the ^Po and the ^Si Di (deuteron) channels. In other 
channels the improvement is moderate, indicating miss- 
ing shorter range contributions to Eq. ([1]). Although a 
deeper understanding on why this dramatic reduction of 
the number of counter-terms happens would very helpful, 
and we have not attempted a large scale fit, it is very re- 
warding that the implementation of the correct and fairly 
complete long range physics deduced from One and Two 
Pion Exchange in conjunction with the requirement of 
renormalizability provides a rather reasonable descrip- 
tion of the NN scattering data below pion production 
threshold. 
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APPENDIX A: ANALYTICAL DETERMINATION 
OF FINITE CUT-OFF CORRECTIONS 

In this appendix we determine the finite cut-off cor- 
rections to renormalized deuteron properties when the 
the two coupled channel potential is so that there is one 
attractive and one repulsive short distance eigenpoten- 
tial. This is the case of OPE and RB-TPE potentials 
discussed in this paper. In this case it is simplest to 
discuss the auxiliary boundary condition 



u(r,) = 



(Al) 



The analysis of other auxiliary boundary conditions such 
as u'{r(.) = which has actually been used in the numer- 
ical calculations is a bit messier but the final conclusion 
is essentially the same as with Eq. (lAip . 

From the superposition principle of boundary condi- 
tions we may write 



u{r) = us{r) +riUD{r) , 
w{r) = ws{r) + -qwDir) , 



(A2) 



where {us{r),ws{r)) and (w_D(r), w^i (r)) solve the 
deuteron problem, Eq. (IM)) . with the long distance 
boundary conditions Eq. ([^ when taking [As.Ad) = 
(1,0) and {AstAd) — (0,1) respectively. On the other 
hand, at short distances the coupled channel potential is 
diagonalized by an orthogonal transformation 



cos 9 sin 9 
— sin 9 cos 6 , 



(A3) 



where the mixing angle 9 depends on the parameters of 
the potential only. For instance, for OPE one has cos 9 = 
— 1/-\/3, see Ref. HH]. The reduced potential behaves as 



U(0 - - G 











G 



(A4) 



where i?+ and i?_ are the Van der Waals length scales 
associated to the repulsive and attractive channels re- 
spectively. Defining the general short distance solutions 
of the decoupled problems 
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2/4 



ai,+ exp 



2 



n/4 



Ci - sm 



n — 2 V r 



2 /i? 



6j,+ exp 



2 /R 



n — 2 \ r 



n — 2 \ r 



•ft 



(A5) 



where i — S, D and 0.;.+ , Ci._ , (/^i are fixed constants 
which depend on the potential only and may be deter- 
mined by integrating in from long distances the asymp- 
totic wave functions {us{r),ws{r)) and {ui:i{r),WDi'r)). 
We remind that n — 7 for the RB-TPE potential while 
n = 3 for the OPE potential. Note that we must include 
here also the diverging exponential at the origin for the 
repulsive eigenpotential. The solutions at short distances 
behave as 



(A6) 



where i — S, D. Thus, using the short distance boundary 
condition, Eq. (jAl|) which selects the regular solution at 
the origin and eventually kills the diverging exponentials 




when Tc 
get 



0, and gathering all subsequent equations we 

usjrc) 
unirc) 

cos9 vs,+ {rc) + sine vs,^{rc) ^^^^ 



cos9 V£)^^{rc) + sin6'w£)._(rc) 



The limiting value is controlled by the short distance di- 
verging exponentials in Eq. (|A5p and is given by 



77(0) = 



bs,+ 
bD.+ 



(A8) 



Deviations from this value for small cut-offs Vc can be 
directly determined from Eq. (jA5[) yielding 



77(0) 



1 + tan ( 



Cd.- 



exp 



2 fR, 

n — 2 \ Vr 



sm 



2 ( rr 



(A9) 



showing that, up to oscillations, finite cut-off corrections 

-1 _i 

in the deuteron are Oie-"^" ) for OPE and ©(e"'^- ) for 
RB-TPE. The generalization to other auxiliary boundary 
conditions and other deuteron properties is straightfor- 
ward with identical result in the order of finite cut-off 
effects. 

The case of scattering states is more tedious and will 
not be discussed in detail here but can also be analyzed 



with a combination of the coupled channel formulas of 
Ref. [i^l (see Sect. V of that paper) and the bound state 
results of the present appendix. Finite cut-off effects for 

the S-matrix scale similarly as in the bound state case, 

_ 1 

i.e. up to oscillations they are C'(e^''^ ) for OPE and 

_ 5 

0(e-''<= ') for RB-TPE. 
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